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ANALYSIS OF STATIC AEROELASTIC BEHAVIOR 
OF LOW- ASPECT-RATIO RECTANGULAR WINGS 


By John M. Hedgepeth and Paul G. Waner, Jr. 


SUMMARY 


Slender-body theory Is used in conjunction with, plate theory to 
analyze the static aeroelastic-divergence behavior of low-aspect-ratio 
rectangular wings of constant thickness when chordwise deformations are 
considered. In the analysis, the spanwise variation of the deflection 
is restricted to a parabola but the chordwise variation is allowed com- 
plete freedom. Results show the variation of the divergence speed and 
mode shape with the aspect ratio. Comparisons are made with results 
obtained by using approximate (linear, parabolic, and cubic) chordwise 
deflection shapes. 


INTRODUCTION 


Methods for predicting the divergence speed of wings when chordwise 
deformations are neglected have been treated extensively in the past. 

For wings with fairly large aspect ratios, accurate results have been 
obtained; however, for wings with low aspect ratios, the chordwise defor- 
mations can no longer be neglected. 

A number of analyses are available that deal with the effects of 
chordwise deformation on divergence . Among them are the works of Miles 
(ref. 1) and Biot (ref. 2). Miles considered the chordwise divergence 
of a delta wing cantilevered along its trailing edge. He assumed that 
the deformations of the wing were cylindrical with straight-line gener- 
ators in the spanwise direction. Biot, on the other hand, considered an 
unswept wing and included both spanwise and chordwise structural effects. 
In Biot's work, the use of aerodynamic strip theory limits the value of 
his analysis for low-aspect-ratio wings. 

In the present paper an analysis is made of the divergence of 
very low-aspect-ratio cantilever plates of uniform thickness (see fig. 1) 
Allowance has been made for the presence of additional discrete chordwise 
stiffe ning elements. Although the analysis includes both spanwise and 
chordwise structural effects, primary emphasis is placed on the chordwise 
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deformations. Indeed, the primary purposes of this study are to deter- 
mine the types of chordwise deflection shapes which can he expected in 
very low-aspect-ratio wings and to assess the accuracy resulting from 
the use of appro xima te chordwise mode shapes in aeroelastic analyses. 

Slender-body theory (see ref. 3> for example) is used to determine 
the aerodynamic loads, and plate theory is used in conjunction with a 
potential-energy approach to determine the deformations . Numerical 
results are presented for wings with various aspect ratios and chordwise 
stiffnesses. 

SYMBOLS 

A aspect ratio 

a-j_, & 2 > a 3 coefficients defined by equations (17) 
c wing chord 

D 
E 

(El) e 

(EI)± 

f(x), f(S) 

N 

Pl(x) 

p(x,y) 

q. 

E coefficient defined by equation (l8) 


plate stiffness in bending, Et^/ 12 (l - p 2 ) 

Young's modulus of elasticity 

total effective bending stiffness of chordwise stiffeners 
(See eq. (19) . ) 

bending stiffness of ith stiffener 

chordwise deflection shape 

total number of chordwise stiffeners 

generalized aerodynamic loading (See eq. (13).) 

lift per unit area, positive in z-direction 

pU 2 

dynamic pressure, — — - 


s 


wing semi span 
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t wing thickness 

U free-stream velocity 

V]_ generalized aerodynamic leading-edge shear (See eq.. (13).) 


V(y) 
w(x,y) 
x,y,z 
a, (3,0 

A 

P 

n 

a 

P 

0(x,y,z) 

Subscript: 


leading-edge shear 

wing deflection, positive in z-direction 

coordinate system (See fig. 1.) 

parameters defined by equations (20) and (23) 

divergence-speed parameter, 

Poisson's ratio 
nondimens ional coordinates 
total potential energy of the system 
dummy variable of integration for y 
free-stream density of fluid 
perturbation velocity potential 


5jt qs^ 
lj-8 D 


x/s, y/s, respectively 


i 


integer denoting stiffener number (See fig. 1. ) 


ANALYSIS 


The wing configuration considered herein consists of a rectangular 
plate of constant thickness with a number of constant-stiffness beams 
in the chordwise direction (fig. 1). The wing is assumed to be supported 
along its center line and only symmetrical deflections are considered. 

The analysis is based on the assumption that the aspect ratio of 
the wing is very low. For such a wing, the deflections will generally 
vary in a much more complicated manner in the chordwise (x) direction 
than in the spanwise (y) direction. It therefore seems reasonable to 
assume a simple spanwise variation for deflection and allow the chordwise 
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variation to be arbitrary. Since the simplest type of spanwise defor- 
mation compatible with the support boundary conditions is a parabola, 
the following assumed deflection shape is used in the analysis: 


w(x,y) = y%(x) (l) 

where f(x) is the chordwise deflection shape. 

The assumption of low aspect ratio also implies that Jones' slender 
body theory (ref. 3), in which streamwise perturbations are considered 
negligible in comparison with perturbations in the crossflow directions, 
can be used to determine the aerodynamic loading caused by the deforma- 
tion of equation (l) . Although slender-body theory is exact only in the 

limit (l - M 2 )a 2 -» 0, it is employed herein because of its eaBe of 
application; the use of slender-body theory leads to a particularly sim- 
ple formulation of the divergence problem and, hence, allows a detailed 
investigation of the divergence behavior of very low-aspect-ratio wings. 

In order to determine the deformations resulting from aerodynamic 
loads, plate theory is employed. The principle of minimum potential 
energy is used to derive the differential equation of equilibrium for 
the function f(x). In applying this principle, the difficulties arising 
from the presence of nonconservative forces are circumvented by treating 
the loads as fixed quantities during the variation. The aerodynamic loadt 
are then substituted into the differential equation and the solution 1 b 
obtained. 


Aerodynamic Forces 

For slender-body theory, the velocity-potential equation for line- 
arized flow reduces to Laplace's equation in the crossflow plane: 


5 2 0 

cty-2 


( 2 ) 


The boundary conditions are 


(x,y,0) 


Uy 2 — 
dx 


(5) 


on the wing and <jb = 0 at the tips. Furthermore, derivatives of the 
potential at Infinity must be zero. 
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It is desired to calculate the potential 0 and hence the pressure 
distribution resulting from the given deformation shape. 

In order for equation ( 2 ) to be satisfied, with the boundary con- 
dition that 0 = 0 at the tips, the following relationship between 

— (x,y,+0) (the limit of 30/dy as z approaches zero from above) 

oy 

and ^ (x,y,0) must exist (ref. h): 
oz 


S0(x,y,+O) 

dy 


it 


■ _ V 2J_ 


8 /s2 _ a 2 50 (x, cr,0) 


s - y 


y - a 


Sz 


da 


00 


Substituting ^ (x,y,0) from equation ( 3 ) and integrating gives 
oz 


0(x,y,+O) = -2 f '(x)^ + y 2 j\/s 2 - y 2 


(5) 


where the prime denotes differentiation with respect to x. 

The lift per unit area of a thin wing in terms of the velocity 
potential is 


p(x,y) = 2pU 

dx 

Thus, the lift per unit area is, by substitution from equation ( 5 ), 

p(x,y) — | qf"(x)^ + y 2 )|/s 2 _ y 2 ( 6 ) 

In addition to the distributed pressures acting on the wing, there 
is a concentrated load acting along the leading edge . This concentrated 
load arises as a result of the application of the slender-body theory to 
the rectangular plan form; that is, the velocity potential is zero ahead 
of the leading edge and then jumps to a finite value at the leading edge. 
The magnitude of the resulting concentrated load is 


V(y) = 2pU0(O,y,+O) 
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which gives, by use of equation (5), 


V(y) 


= # ’(0) 
5 




(7) 


Structural Equilibrium 


The total potential energy of the system under consideration con- 
sists of the strain energy of the wing plus the potential energy of the 
aerodynamic forces. The strain energy of the wing is made up of the 
strain energy of the plate (ref. 5) 




2 VL &k<&L + 

^ dx2 £y2 + 


13 + 2(1 



dy dx 


and the energy of bending of the chordwise stiffeners (N in number) 


if (EI,i / 0 ' 

i=l 



2 

dx 


In these expressions, D is the plate stiffness (^Ety 12(1 - p 2 )) , (El) 
is the bending stiffness of the ith stiffener, and yj. is the spanwise 

coordinate of the ith stiffener. (Note that because of symmetry, only 
half of the wing is considered.) 

The foregoing analysis of the aerodynamic forces has indicated that 
the wing loads are composed of two parts, a distributed lift p(x,y) 
and a concentrated load V(y) located at the leading edge. The poten- 
tial energy of the lift per unit area is 



p(x,y)w(x,y)dy dx 


and the potential energy of the concentrated load is 
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The total potential energy of the system is the sum of these dif ferent 
energies, or 


c o s 


11 ■ i/o So 




dy dx + 


(EI)i /q C ^ “ IJi: p(^y) w Cx,y)dy - 


i=l 
, s 


f V(y)w(o,y)dy 

Jo 


( 8 ) 


Substitution of the expression for w(x,y ) given "by equation (l) 
into equation (8) results in the following expression for the total 
potential energy: 


H 


- 8 /„ 


|5 + i 

5 D 


N 


(El)jyi 


i=l 


f" 2 + W 2 + ^ s 5 uff" + 
3 


| s3(l _ n)f >2 
3 


dx - T C P 1 (x)f dx - V-^ 

J 0 


f(0) 


where 


(9) 


r 3 

Pl(x) = / p(x,y)y2dy 

■JO 


Vi 


r 3 

= / V(y)y‘ 

Jo 


'dy 


(10) 
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Minimization of the potential-energy expression hy means of the 
calculus of variations gives 


&JI = 0 



(EDiMi 4 


Sf" 2 + ls8f 2 



s3n8(ff") 


+ 


8 

3 


b3(i - 


|i)5f 




P 1 (x)Sf dx - V^ffO) 


Integrating by parts and collecting terms resultB in the differential 
equation 


s£ 

5 


+ 



(EI) 1 y 1 1< - ■ § ! s 3f" + 4sf 


Pl(x) 

D 


( 11 ) 


and the following boimdary conditions: 
At both x = 0 and x = c, 


T + S X (EI,iyi “ 

1=1 


f" + I s 5 pf 

3 


= 0 


At x = 0, 


N 

T + D X (EI) i y i 4 


i=l 


f'" _ (8 - lOp) s 3f< 

3 


^-0 


(12a) 


(12b) 
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At x = c, 


f + 1 X (Ei,iyi1 

1=1 


f m - £ 8 “ 10 K) s 3 f > = 0 
3 


(12c) 


Aeroelastlc Solution 


The differential equation and "boundary conditions which describe 
the deflection of the structure due to the generalized aerodynamic loads 
P^(x) and Vj_ are given by equations (ll) and (12) . The values of the 

generalized loads can be obtained by substituting p(x,y) and V(y) 
(eqs. (6) and ( 7 )) into equations (10). This substitution gives 


P l( x ) = “i2 s6 1 f "( x ) 
Vi s 6 qf'(0) 


1 




(13) 


If these values of Pj_(x) and V-j_ are substituted into the differential 

equation (ll) and the four boundary conditions (eqs. (12) ) , and the equa- 
tions are nondimensionalized by letting 


i 



and 


Hi 


s 


the differential equation becomes 


fW (l) + |(A - a^f'd) + ^ f(|) = 0 
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The four boundary conditions became: 
At |=0 or 


f"(|) + f(|) = 0 


At |=0, 


At 



where 


and 


f m (i) + 


4 A 

\ a 2 

LhC 

■ a v “ 


f ' ( 6 ) = 0 


f"U) - | a 5 f'(|) = 0 



( 15 ) 


(16) 


( 17 ) 


5(EDe 


R = 1 + 


sD 


(18) 
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In equation ( 18 ), (El) e is the total effective bending stiffness of 
the chordwise stiffeners and is defined as 

N 

(El) e = (EI )^ 4 (19) 

i=l 


The characteristic roots of the differential equation are 


- a i) ± ]/(* - a i ) 2 

The roots can be seen to be dependent on the value of - .a^j . For 

- a ^ 2 > 5 R, all the roots are imaginary^ for (\ - a ^ 2 < 5 R, there 

are two pairs of conjugate complex roots. Consider first ^A - a j ^ 2 > 5R» 

For this case, the characteristic roots of the differential equation can 
be written 


±i(<x ± p) 


where 


0 - l/i ||( x - a i) + t® 

13 ' l/F |/( x - a i) - ® 


( 20 ) 


The solution to the differential equation is then 


f(|) = C-l cos(a + p)f= + C 2 sin(a + 0 )g + 
C 3 cos(a - p)| + sin(a - 0 )§ 


( 21 ) 
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The value of f(g) is now substituted into the boundary conditions 
given by equations ( 15 ) ; this substitution yields a set of homogeneous 
equations in the constants C^, Cg, C^, and C^. The existence of a 

nontrivial solution requires that the determinant of the coefficients 
be equal to zero, or 


b^ cos(a + bj_ sin(a + p)j| 1*2 cos(a - p)£ b£ sin(a - 

0 bj 0 b^ 

-b 5 sin (a + p)-| b^ cos(a + p)-| -bg sin(a - p)| bg cos(a - p)§ 


where 


= 0 


*1 = ~ - (a + P Y 


R 


b 2 - ^ - (a - p ) 2 


= (a + P) 

bjj. = (a - p) 

b^ = -(a + p) 
bg = -(a - p) 


(a + p) 2 + ^ - |(A - aj 


(a - P ) 2 + ^ - i(X - aj) 


^ a j + (a + p)^ 


! ] 

^ a^ + (a - p) 2 J 


After the determinant has been expanded and simplified, the following 
equation is obtained: 


sin^p — 
r s 

/ \2 
h) 

(a - a-L + A-l)(a - a-|_ + Ag) 

sin2a ■§ 

D 

\ a / 

(“^ + a i + a 3)(^ - a l “ A k) 


( 22 ) 
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where 


Al ' 2 VI ( 5 " w) " 2 a 2 

a 2 = 1 \^(i + is) + ^ - f|) 

A 3-|/f( 5 -i-) + l a 2 

^ - | ^ + S^) - ( a 5 - I) 

Equation (22) is valid only when ^A - a-j^ 2 = 5R. For values of 

(a - aq) 2 < 5R, the quantity p becomes imaginary. The corresponding 
divergence equation_can be obtained directly from equation (22) by 
replacing p by ip where 


5 - #1/^ - (* - -j) 


(23) 


The result is 


sinh 2 j3f /^2 

(A - a-^ + AjJ^A - a-^ + 

sin 2 a § ^ 

^-A + aq + A'j^A - aq - Alj.^ 


(2k) 


Equation (22) (or (2*0) gives an implicit relationship between the 
aerodynamic parameter A and the structural parameters (El) e ysD, \±, 

and c/s. If (El) e /sD and n are fixed, then the values of A n are 
known and the variation of A with c/s can be computed by a trial - 
and-error process. Since c/s appears in equation (22) (or ( 2 ^)) more 
simply than A, this computation can be accomplished most easily by 
choosing values of A and solving for the corresponding values of c/s, 
where they exist. 
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Approximate Solution 

In the aeroelastic analysis of low-aspect-ratio wings of practical 
construction, a closed-form solution such as that described in the pre- 
ceding section would not be feasible; some sort of further approximation 
would be necessary. One type of approximate procedure is to restrict 
the allowable chordwise deflection shape to a polynomial of finite degree. 
The present conf iguration (fig. l) furnishes an excellent test for the 
accuracy of this procedure . 

The consequences of assuming linear chordwise deflections can be 

( 0 &) \ 

obtained by allowing the chordwise stiffeners to be rigid l = . 

\ sD / 

In order to obtain results for higher degrees, it is necessary to start 
with the energy expression (eq. (9))* Letting 

f(*> - a 0 ♦ *i(f) ♦ a 2 (ff ♦ a 3 (f) 3 ♦ . . . 

then minimizing the potential, energy with respect to d-j_, dg, d^ . . ., 
and finally substituting the appropriate expressions for and Vj_ 

from equations (13) yields a set of homogeneous simultaneous equations. 
For the cubic approximation, the equations are: 


5 

1(2.) + A. 

2W + c/s 

Bg + 2 A 

b 5 + 

3A - 
s 

^0 


0 

5/2 

Bl 

Bj + A 

B5 + 

2A - 
s 

di 


0 

B 2 

«L. 

0 

Bj, + =. A 

-f- 


d 2 


0 

(c/s ) 2 

c/s 


( 

2 s 



b 5 

B 6 

+ B A 

Bq + 

|a J 



0 

(c/s ) 5 

(c/s ) 2 

c/s 2 

5 5 

j 




( 25 ) 
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where 



The value of A can now he found for known values of (EI) e y^sD, p., 

and c/s by setting the determinant of equation (25) equal to zero. 

For the parabolic approximation, the corresponding determinant can be 
found by deleting the last row and column of the matrix in equation (25). 

Computed results as obtained from this approximate analy sis as well 
as those obtained from the more exact analysis are presented and dis- 
cussed in the following section. 



16 


NACA TN 3958 


RESULTS AND DISCUSSION 


The results of the exact divergence-speed calculations for a low- 
aspect-ratio wing for which Poisson's ratio equals 1/3 are shown in 
figure 2. The terms "exact" and "approximate" are used to denote exact 
and approximate solutions to the approximate formulation of the static 
aeroelastic problem presented herein. In this figure (fig. 2) the 
results are given in the form of a plot of the divergence- speed param- 
eter A against c/s ’for various values of the beam-to-plate stiffness 
ratio (El) e /sD. Results are shown only for values of c/s greater 

than unity because of the obvious limitations of the slender -body theory 
for smaller values. As is to be expected the effect of adding chordwise 
stiffness to the wing is to increase the divergence speed. The values 
of A for (El) e ^sD = 00 (corresponding to a wing in which the chord- 

wise bending is neglected) increase without limit as c/s becomes large 
(or aspect ratio becomes small). On the other hand, the values of A 
for a finite (El) e /sD apparently reach a constant value. The curves 

actually undulate slightly; however, the deviation from a straight line 
is so small that it is not apparent in the figure. 

The accuracy of the linear, parabolic, and cubic approximations 
for the chordwise deflection shape is illustrated in figure 3. In this 
figure the variation of A with c/s for the approximate solutions is 
compared with that for the "exact" solution for the case of zero chord- 
wise stiffening. It can be seen that the linear approximation is unsat- 
isfactory for all values of c/s greater than unity. The parabolic 
approximation gives values of A accurate to within 10 percent of the 
exact value for c/s less than 2.2. The range of this accuracy is 
extended for the cubic approximation to c/s equal to 4. 

An examination of the mode shapes at divergence is of Interest. 

The mode shapes for (EI^^sD = 0 are shown in figure 4 for three 

values of c/s. These values of c/s correspond to the positions indi- 
cated by the ticks on the curve for (EI) e /sD =0 in figure 2. It can 
be seen that the effect of increasing c/s is essentially to extend 
the mode shape rearward; the shape near the leading edge changes only 
slightly. 

Also shown in figure 4 is the mode shape for c/s = ». This shape 
can be obtained in the following manner: 

For all the solutions obtained for finite c/s, the parameters are 
such that the value of 0 (eq. (20)) is imaginary; therefore, two of 
the homogeneous solutions of the differential equation (l4) approach 
zero and two approach infinity as £ approaches infinity. If the latter 
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two solutions are omitted, an equation governing divergence is obtained 
through the use of the boundary conditions at the leading edge only. 
This equation is 


A = a-^ + A j ( 26 ) 

which gives an asymptotic value that agrees within plottable accuracy 
with the flat portions of the curves in figure 2. The corresponding 
mode shape is 


f(e\ e“^sln(o4 + 6) 
sin 0 


(27) 


where 


0 = tan 


-1 


2o£ 


is . a 2 + 52 

R 


Note that equations ( 26 ) and (27) hold for finite values of (El) e /sD 
as well as for (El) e /sD = 0. The variation of A for large c/s with 
(El) e jsh is shown in figure 5 . The increase in divergence speed resulting 
from chordwise stiffeners is clearly evident. 

CONCLUDING REMARKS 


The present analysis of the static aeroelastic divergence of low- 
aspect-ratio rectangular wings indicates that the deflection shape in 
the chordwise direction has an increasing number of waves as the aspect 
ratio is reduced. The inclusion of additional chordwise stiffening 
reduces the amount of chordwise bending and, consequently, increases the 
divergence speed. Approximating the chordwise deflection shape by para- 
bolic and cubic curves yields divergence speeds in fair agreement with 
the predictions of the more exact theory if the aspect ratio is not too 
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low. The cubic approximation Is more accurate than the parabolic one 
primarily in that configurations of lower aspect ratio can be treated 
without the error becoming excessive. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., January 24 , 1957. 
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Figure 2.- Divergence -speed results obtained by "exact" analysis for 

,=1/3. A. 12 i|(|) 5 . 




Figure 3 .- Comparison of results obtained with approximate chordwise 
deflection shapes with "exact" results for (El) e /sD = 0. 







